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ABSTRACT 



We study the thermal formation history of four simulated galaxies that were shown in iNaab et al.l 
(|2007t l to reproduce a number of observed properties of elliptical galaxies. The temperature of the 
gas in the galaxies is steadily increasing with decreasing redshift, although much of the gas has a 
cooling time shorter than the Hubble time. The gas is being heated and kept hot by gravitational 
heating processes through the release of potential energy from infalling stellar clumps. The energy is 
dissipated in supersonic collisions of infalling gas lumps with the ambient gas and through the dynam- 
ical capturing of satellite systems causing gravitational wakes that transfer energy to the surrounding 
gas. Furthermore dynamical friction from the infalling clumps pushes out dark matter, lowering the 
central dark matter density by up to a factor of two from z — i to z — Q. In galaxies in which 
the late formation history {z < 2) is dominated by minor merging and accretion the energy released 
(E' ~ 5 X lO^^ergs) from gravitational feedback is sufficient to form red and dead elliptical galaxies by 
z ~ 1 even in the absence of supernova and AGN feedback. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: formation — galaxies: evolution 
— methods: numerical 



1. INTRODUCTION 



In the standard cold dark matter (GDM) picture of 
galaxy formation, gas falling into dark matter halos 
is shock-heated approximately to the halo virial tem- 



perature, Tvi 



10" 



c/167kms ^) K maintain- 



ing quasi-hydrost atic equilibrium with t he dark mat 
ter comp o nent. (iRees fc Ostrikeij Il977t iBinnevI [19771 : 
Isiigfl977t iWhite fc Reed 1X9781) . The gas will cool, los- 
ing its pressure support, and settling into a centrifu- 
gally supported di sk while conserving its specific an- 
gular momentum (Fall & Efstathiou 19 801). However, 
recently Birnboim fc Dekel (2003 ): Dekel fc Birnboiml 
(|2006( l showed using one-dimensional models that galaxy 
halos can only shock-heat infalling gas if the cooling rate 
for gas behind the shock is lower than the compression 
rate of the infalling gas. This criterion translates to 
a roughly redshift independent critical minimum mass 
for halos that are able to shock-heat the infalling gas of 
-^^shock ~ 10^^ '^M©. Less massive halos are not able to 
support stable shocks and the majority of their gas is ac- 
creted cold. However, in the general three-dimensional 
case galaxies above the critical shock mass can also be 
fed with cold gas along filaments penetrating deep inside 
the hot halo (jOcvirk et al.ll2008f ). 

Recent simulations also predict that the baryonic 
growth of most galaxies is dominated by smooth accre- 
tion events of either hot or cold gas, as opposed to in- 
divid ual merger events (|Genel et a]||2008t iBrooks et al.l 
120081 ). Observations of massive star-forming galaxies at 
z ~ 2 have found a preponderance for thick gas-rich ro- 
tating disks with large tur bulent motions and only a mi- 
nority of mai or me rgers (jForster Schreiber et all l2006t 
IShapiro et al.ll2008[ ). The observed high-redshift disks 
exhibit large turbulent motions that could potentially 
be the result of gravitational ener gy release from col d 
accretion flows feeding the galaxy ( Genzel et al.|[2008[ ). 



Given the bimodal temperature distribution of the ac- 
creted gas and the characteristic mass scale of A/ghock, it 
is tempting to link this to the observed bimodality in the 
galaxy population dividing them into a red and blue se- 
quence with a critical stellar mass of Merit — 3 x IO^^Mq, 
which was reve aled in recent large gal axy surveys (e.g. 
iBell et al.l[2003HKauffmann et al.ll2003l ). Galaxies below 
this critical mass are typically blue, star-forming disk 
galaxies, whereas galaxies above Merit are dominated by 
red spheroidal systems with old stellar populations. 

The observed bimodality can be explained if star 
formation in halos above a critical threshold mass of 
M ~ W^Mp, \s su ppressed (e.g. iBower et all 120061: 
ICattaneo et al.l[2b06f l. The quenching mechanism needs 
to be both energetic enough to trigger the quench- 
ing and long-lasting enough to maintain the quench- 
ing over a Hubble time. In addition to the quench- 
ing by sh ock-heated gas abov e a critical ha lo mass 
(jDekel fc Bir nboim 2006 : Birn boim et al.ll2007f ). poten- 
tial qu enching mechanisms inc lude the feedback from 
AGNs (jCiotti fc Ostrikeii [200l . gaseous major merg- 
ers triggering st a r burst and/or quasar activity (e.g 
Naab et all 120061: iHopkins et all l2007l : [Johansson et al.l 
20091) aiidgravitational quenching by clumpy accretion 



( Dekel fc Birnboi"mll2008f l. Environmental effects and the 
input of gravitational energy from infalling clumps are 
real and automatically included in high-resolution nu- 
merical simulations. In semi-analytic models these envi- 
ronmental effects are not u sually included (see, however, 
iKhochfar fc Ostrikeill2008l for an exception.) 

In this letter we study in detail the thermal forma- 
tion history of a set of very high resolution cosmological 
re-si mulations of indiv idual galaxies. In an earlier pa- 
per (jNaa b et al.l [2007!, hereafter N07) we showed that 
the simulated galaxies reproduce a number of observed 
properties of elliptical galaxies, although the simulations 
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Fig. 1. — The evolution of the temperature profile (left panel), the entropy distribution (middle panel) and the phase-space diagram 
(right panel) of all gas within rvir as a function of redshift (top to bottom) for halos Al (black), C (green), E (orange) and halo A2 (red), 
where /jiot is the fraction of gas with T > 2.5 X 10^ K. Even though much of the gas has a cooling time shorter than the Hubble time its 
temperature steadily increases due to the enumerated gravitational heating processes. 

include neither stellar nor AGN feedback. We show that 
the gas is heated and kept hot by gravitational feed- 
back, efficiently terminating star formation and leading 
to the formation of systems with evolved stellar popula- 
tions and low star formation rates by ^ ~ 1. Gravita- 
tional feedback comes in many forms. Incoming lumps 
or streams of cold gas ultimately come to rest deposit- 
ing their potential energy frictionally (i.e. through de- 
caying turbulence). Supersonic collisions of infalling gas 
with the ambien t gas lead to propagating shock waves 
(|Rvu et al.ll2003[ ) which also deposit entropy throughout 
the system. Gradual accretion to the gaseous envelope 
of the galaxy adds weight, causing contraction and PdV 
work on the ambient gas. And finally infalling satel- 
lite systems captured through dynamical friction cause 
gaseous wakes fr om which energy is transferred to the 
surrounding gas (jOstrikeij[r999 ). These are all variants 
of the processes by which gravitational energy released 
by infalling matter can be transferred to the gas main- 
tained in quasi-hydrostatic equilibrium within galaxies, 
heating it and tending to counteract radiative losses. 



2. SIMULATIONS 

The ACDM initial conditions assu me scale-invariant 
adiabatic fluctuations with T — 0.2 (jEfstathiou et al.l 



119921). Throughou t this letter we use a WMAP-1 
(jSpergel et al.l l2003f ) cosmology with a slightly lower 
Hubble parameter oi h — 0.65^ with CT8=0.86, /& — 
^^h/f^m=0.2, f2o=0.3, and Ao=0.7. The galaxies were 
simulated at high resolution usin^ the volume renormal- 
ization technique ()Katz fc White" iggS') by selecting tar- 
get halos at 2; = in low-density environments from 
a low-resolution dark matter simulation. In the re- 
simulations the particle number of the gas and dark mat- 
ter particles was increased to 100^ (Halos Al, C and E) 
and 200'^ (Halo A2) within a cubic volume at redshift 
z — 2A containing all the particles that ended up within 
the virialized region of the halos at z = 0. The virial 
properties of the galaxies at z = are summarized in 
Table [H 

The simulations were performed using the TreeSPH- 
code GADGET-2 (Springel 2005) on shared-memory ma- 
chines hosted in Cambridge, Munich and Princeton. The 
code includes star formation and optically thin radia- 
tive cooling for a primordial c omposition of hydrogen 
and helium ( Katz et al.l 119961 ). We included a spa- 
tially uniform redshift-dep endent UV background radia- 
tion field with a modified (|Haardt fc Madaulll996[ ) spec- 

1 h defined such that Ho=imh kms-^Mpc-l. 
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TABLE 1 
Galaxy properties within rvir 



Redshift 

0.1 0.3 0.6 1 1.5 2 346 



Galaxy 




-^-^stai's 




M-OM 


r ■ 

' vir 




Al 


169 


21.7 


13.3 


134 


258 


270 


CI 


125 


17.2 


9.3 


98.5 


233 


258 


El 


134 


19.0 


9.3 


106 


239 


307 


A2 


179 


22.8 


11.6 


145 


263 


232 



"Total masses M in lO^'A/©. 

'Virial radius in kpc, defined as the radius enclosing an overdcn- 
sity of 200 times the critical density Pcrit- 
'^Maximum circular velocity in km/s. 

trum. We modify th e self-r egulated feedback model of 
ISpringel &: Hernguistl (|2003f) by turning off the feedback 
from supernovae, eliminating the two-phase description 
of the ISM of star-form ing particles with p > Pthres h 
(jMcKee fc Ostrikedll977HJohansson fc Efstathimjll2006f ). 
and thus transforming their entire gas component into 
the cold phase with a pre-set temperature of Tcoid = 
1000 K. The star formation time scale is set by = 
*o(Pthresh/p)^^^, whcrc we set the free parameters to 

= 1.5/i"iGyr and pthrosh = 7 x lO'^^h'^ g cm^^. 

The gravitational softening length for the 100^ (200^) 
runs was fixed in comoving units of egas = e^, = 
0.25 (0.125) kpc and edm = 0.5 (0.25) kpc until z ^ 9, 
after which the softening remained fixed in physical coor- 
dinates at the same values. Our lowest (highest) resolu- 
tion runs have particle masses of mgtars = 1.05 x 10^ Mq 



(mstars = 1.3 X IO^Mq), mgas = 2.1 X IO^Mq (mgas = 

2.6 X IO^Mq) and mdm = 8.2 x 10*^ Mq (mdm = 1-03 x 
1O^M0) for the stars, gas and dark matter, respectively. 



3. GRAVITATIONAL HEATING 
3.1. Temperature structure of the gas 

We summarize the evolution of the gas temperature 
at redshifts z = 0, 1,4 for our model galaxies in Fig. [T] 
by showing the temperature profiles (left panel) , the en- 
tropy distributions (middle panel) and the phase-space 
diagrams (right panel) of all gas within the virial radius. 
The temperature of the gas is increasing in all galaxies 
with decreasing redshift reaching T^lO^Katz^O, 
with cool, star-forming gas only found at the very centers 
of the galaxies. In the middle panel we plot the distri- 
bution of the gas entropy, defined as 5 = kTn'"^/^ and 
define the cooling time of the gas as 



S 1.5(/ie/M)' ■ (10 keVcm2)3/2 



10 keVcm2 



{kT)^/^A{T, Z) 



(1) 



where k is the Boltzmann constant, A(T, Z) is the cool- 
ing function, and /i ~ 0.59, — 1-1 for a fully ionized 
gas. The first term is a measure of the entropy and there- 
fore an adiabatic invariant, whereas the second term only 
depends on the temperature (T) and metalicity (Z) o f 
the gas (|Scannapieco fc Oh|[200llKhalatvan et al.ll2008D . 
The second factor in Eq. [T] has an absolute minimum, 
corresponding to tcooi ~ 2 Gyr for our primordial cooling 
function (and S=10 kevcm^, with different 5* correspond- 
ing to different icooi)- We plot in Fig. [T] the entropy 
values corresponding to minimum cooling times of 0.1, 
1 and 10 Gyr as dashed lines. At high redshifts the en- 
tropy distribution of the gas is bimodal with cold, high- 
density, star-forming gas forming a low entropy peak and 
lower density, hot shock-heated gas forming a high en- 
tropy peak. At lower redshifts the available low entropy 
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Fig. 2. — The cumulative change in binding energy AiJbind 
summed from z=5 to z=0 separately for in situ (circles) and ac- 
creted stars (squares) as a function of lookback time (redshift). 

gas has been consumed by star formation and we are 
primarily left with dilute shock-heated gas with cooling 
times of the order of tcooi ~ 5 — 10 Gyr. 

At high redshifts (z > 4) the gas flows into the galaxies 
predominantly in the cold phase with the hot gas (defined 
as T > 2.5 X 10^ K) fraction being low at /hot < 30%. 
A large fraction of the gas can also be found in the cold 
star-forming phase. With decreasing redshift the frac- 
tion of hot gas steadily increases; by z ~ 3 the hot gas 
fraction is around /hot ~ 50%. The transition from a 
cold gas accretion mode to a hot gas dominated accre- 
tion mode occurs at z ^ 2 — 3, corresponding to halo 
masses of Mhaio = 3 — 5 x IQ^^Mp^ in g ood agreement 
with the prediction s of lDekel fc BirnboimI (2006). sec also 
iKeres et alj ()20Q5[ ). Below these redshifts the galaxy ha- 
los are massive enough to support stable shocks and most 
of the accreted gas is shock-heated close to the virial 
temperature of the halos. As this transition is mass- 
dependent, it occurs somewhat earlier in the more mas- 
sive halo Al, compared to the slightly lower mass halos C 
and E. At z = the hot gas fraction is /hot ^ 97%, with 
only some residual cold star-forming gas remaining, as 
most of the original gas has either formed stars or been 
shock-heated to high temperatures. 

3.2. Gravitational energy input 

The diffuse gas is heated by gravitational feedback in 
the form of gravitational energy release from infalling 
stellar and gaseous clumps. We quantify this effect in 
Fig. [5] which shows the cumulative change of the bind- 
ing energies (S(Ai?bind)) of in situ (circles) and accreted 
(squares) stars within a fixed radius of r < 30 kpc. 
In situ stars are defined as stars that were born from 
gas in the galaxy (r < 30 kpc), whereas accreted stars 
have formed outside the galaxy and have been accreted 
later on. The I]ins(Ai?bind) reflects the overall mass 
growth and deepening of the potential well, whereas 
I]acc(A£'bind) dcscribcs the release of gravitational en- 
ergy from accreted structures. All of the energy released 
by the growth of the potential energy of the accreted 
stars (i?bind ~ 5 x lO^^ergs) goes into heating the gas 
and moving outwards dark matter and previous gener- 
ations of stars (on average ~ 5 — 10% of the released 
energy ends up as kinetic energy of the stars), whereas 
much of the energy release associated with the potential 
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Fig. 3. — The net heating (solid line) and net cooling rates (dashed lines) for non-starforming (p < Pthrcsh) g^is within rvir at redshifts 
0, 1, 3 and 5 for the four halos. The fraction /gf of dense starforming gas (p > pthreah) is also given. The heating rate dominates over the 
cooling rate at all redshifts for the low-density non-starforming gas. 



energy growth of the in situ stars is radiated away. 

Ah galaxies assemble rapidly at high redshifts (z — 
5 — 3) with S]A(i?bind) of both the in situ and accreted 
stars increasing by an order of magnitude. The change 
in the accreted EaccC^-E-bind) mirrors the mass accretion 
history of the galaxies with the 3:1 merger for Galaxy C 
at z=0.6 and the nearly 1:1 merger for Galaxy E at z=1.5 
clearly visible (see Fig 7. in N07 for details). At 2; < 1 
the integrated change in the potential energy is domi- 
nated by accreted stars for halos A1,C,A2 (dissipation- 
less formation) and in situ stars for halo E (dissipational 
formation), with the binding energy/atom increasing by 
~ 20% from z = \ \.o z = O.The Eacc(Ai;bind) is 2-3 
larger compared to Sins(A£^bind) for halos Al,C and A2 
at z < 1, whereas halo E is dominated by in situ star for- 
mation and I]i,is(Ai?bind) is a factor of four larger than 
Eacc(Ai?bind)- Wc cxpcct gravitational feedback effects 
to be important in the galaxies assembled through dis- 
sipationless processes (Al,C,A2), whereas gravitational 
feedback plays a minor role in the dissipational forma- 
tion of galaxy E. We also see a trend with increasing 
resolution, the final Eacc(Ai?bind) of halo A2 is higher 
by 50%, whereas the Eins(Ai?bind) is lower by 30% com- 
pared to halo Al. Thus we expect stronger gravitational 
feedback effects with increasing resolution as we are able 
to resolve smaller structures. 

3.3. Heating of the gas component 



The energy input from gravitational feedback scales 
with the mass fraction of each component, with roughly 
ftb/i^rn + ^b) = 1/6 and f7„/(17,„ -I- flf,) = 5/6 going 
into the baryonic and DM components, respectively. Ac- 
creted lumps and cold streams of gas deposit energy into 
the gas through dissipation of turbulence. The infalling 
clumps have typical velocities of Vz=o ~ 400 kms~^, 
Vz=3 ~ 200 kms~^, Vz=5 ^ 100 kms~^ as estimated from 
the difference of the escape velocity at rvir and O.lrvir- 
The corresponding sound speeds of the ambient gas are 
Cz=o ^ 150 kms~^, Cz=3 ~ 100 kms~^, Cz=5 ~ 20 kms~^ 
resulting typically in weak shocks with Mach numbers of 
2-5. 

We quantify this effect by studying the terms of the 
entroDV equation of G ADGET- 2 for a given particle i 
(|Springel fc Hernquistl[200l ). 

^ = -^A(p.,uO + il^^m,n.,Vij . ViWij. (2) 
Pi 2 pI ^ 

Here the first term depicts the external radiative cool- 
ing (or heating) of the gas and the second term gives the 
generation of entropy by artificial viscosity in shocks and 
where the internal energy is defined asu = A/{^—l)p^~^. 
In Fig. [3] we plot the cooling and heating rates derived 
from Eq. [2] for all non-starforming gas (p < pthrcsh) for 
our four galaxies at redshifts z=0,l,3,5. The fraction of 
dense starforming gas (p > ptiucsh) at T = 1000 K that 
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Fig. 4. — Cumulative dark matter mass distribution within r < 
2 kpc as a function of lookback time (redshift). The circles give 
the DM mass in our full baryons+DM (BDM) runs and the squares 
the DM mass in a pure DM (PDM) simulation of halo A2. 

is not participating in the cooling/heating process is also 
given in the figure. The shock-induced heating rate of 
the diffuse gas is larger at all redshifts compared to the 
cooling rate. At large redshifts a substantial fraction of 
the gas ('^ 2/3 at z = 5) and (^^ 1/3 at z = 3) is star- 
forming, whereas at z < 1 this fraction is below a few 
percent. At z < 1 the energy input from gravitational 
feedback is able to keep the gas hot in the galaxies form- 
ing dissipationlessly (A1,C,A2), whereas the dissipation- 
ally forming galaxy E exhibits much lower heating rates 
at z = 1. 

3.4. Heating of the dark matter component 

In addition to heating the gas the gravitational feed- 
back heats the dark matter (DM) of the galaxies caus- 
ing it to expand outwards. We quantify this effect in 
Fig. m where we plot the central DM mass within 
r < 2 kpc in our four galaxies (circles) as a function 
of lookback time. The employed gravitational softening 
(e = 0.25 — 0.5 kpc) is a Plummer equivalent Spline soft- 
ening length, for which the force is exactly Newtonian 
beyond h = 2.8e. Thus r < 2 kpc is resolved by 4 — 8 x e, 
which should be sufficient to follow the DM mass evolu- 
tion reliably. 

In the galaxies that form dissipationlessly we see an 
initial increase in the central DM due to adiabatic con- 
traction. After a peak DM mass at z = 3 the central 
DM mass steadily declines due to heating from infalling 
clumps, with the final central DM mass being a factor 
of 2 lower for galaxies Al,A2 and 50% lower for galaxy 
C. In contrast galaxy E shows a constant evolution of 
the central DM mass, with the z—0 and z—3 central DM 
masses being roughly the same. We also plot for com- 
parison the central DM of a pure DM only simulation of 
halo A2 at 200^ resolution (red squares). For the PDM 
simulation we see neither the adiabatic contraction phase 
nor the pushing out of DM, thus resulting in a constant 
central DM mass throughout the simulation. 

This has important implications for the recent esti- 
mates of central dark matter c ontent of elliptical galax- 
ies (e.g. iThomas et al.l [2009t) . The gravitational en- 
ergy release from infalling clumps might help in trans- 
forming an initially cuspy DM profi le to a cored profile 
as seen in our simulation s (see also lEl-Zant et al.l [200]] : 
iRomano-Diaz et al]|2008l l2009f ). Furthermore, once a 



cored DM profile with a hot low-density gas halo is in 
place, it can be maintained through the infall of cold gas 
clouds e ven without the inclusio n of additional feedback 
sources (jKaufmann et al.]|2008f) . Finally, the accretion 
of minor mergers can potentially explain the strong size 
evolution observed in elliptical galaxie s between z = 2 — 
(jvan Dokkumll2008HNaab et al.ll2009f) . 



4. DISCUSSION 

In this letter we have studied the formation history of 
four galaxies that reproduce a number of the observed 
properties for elliptical galaxies. We have shown that 
for three galaxies (Al,C,A2) for which the late forma- 
tion (z < 2) is dominated by dissipationless minor merg- 
ing that the energy release from gravitational feedback 
alone is sufficient to make these galaxies red and dead by 
z = 1, the fourth galaxy (E) has a dissipational forma- 
tion history with significant ongoing star formation at 
low redshifts. Stellar and AGN feedback are both real 
and important, especially in driving winds and outfiows 
that enrich the intergalact ic medium with metals (e.g. 
lOppenheimer fc Davell2006f ) and potentially in determin- 
ing the lower and uppe r cutoffs in the observed galaxy 
luminosity function (B enson et al.l [20031 ). Furthermore, 
the baryon/DM fraction in our simulated galaxies is a 
factor of two too high for their mass compared to pr edic- 
tions from recent lensing data (jMandelbaum et al.l l2b06). 
However, in order to isolate the effects of gravitational 
heating from infalling clumps, which would occur even if 
there were neither stellar nor AGN feedback, we chose to 
run our simulations without any additional energy feed- 
back. 

The efficiency of gravitational feedback can be esti- 
mated using the following simple argument. 



d4 

dr 



rt3 



r r 



1 / ' vir s. 
:l0g( — ) 



(3) 



where the final change in potential energy results from 
integrating from Tvir to a fiducial infall radius r. Gravita- 
tional feedback energy (Ai?)grav — Am*A0(r) resulting 
in 



(A£;)g,av = Am*<log(100) 



4.5 X 10"%3oom*c^ (4) 



where r^-^/r = 100 and wsoo is the circular velocity in 
units of 300 km/s. Inserting the values for galaxy Al at 
z = (u = 271 km/s and to* = 1.23 x lO^M©) into Eq. 
m results in (AiiJ)grav ~ 8 x 10^^ ergs which is remarkably 
similar to the values that can be read off Fig. [2l More im- 
portantly, this energy is of the same order as the feedback 
energy from supernova II (AiJ)sNii ~ 2.8 x 10~^to*c^ 
(for a Salpeter IMF and 10^^ ergs per supernova) and 
AGN feedback energy (Aii;)AGN 5 x IO^'^to^c^ (for 
msH/'m* = 10~^ and a total AGN feedback efficiency of 
0.5%). However, unlike SNII and AGN feedback energy, 
gravitational feedback energy scales with (Ai?)grav oc v"^, 
making it proportionally more important in massive sys- 
tems with large circular velocities. Furthermore, the en- 
ergy released from gravitational feedback is distributed 
throughout the galaxy halo, whereas the feedback energy 
from SNII and AGN feedback is confined to the central 
region of the halo in regions with high density cold gas 
that can efficiently radiate away the energy. 
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The process presented here is general, but in detail 
might depend on the implementation of additional feed- 
back. The inclusion of supernova feedback might be im- 
portant in lowering the number of stellar clumps thus 
decreasing the amount of gravitational feedback in the 
system. Another important factor not included in our 
simulations is metal cooling that is expected to increase 
the cooling rates by a factor of ten or more. However, the 
increased cooling rates might to some extent be counter- 
acted by efficient stellar and AGN feedback, especially in 
the central regions of the galaxies. Finally, the derived 



heating rates at the peak of the gravitational feedback 
phase (z = 1 in Fig. [3]) are typically 15-20 times higher 
than the cooling rate of gas with primordial composi- 
tion and thus gravitational feedback energy is expected 
to be an important energy source in the dissipationless 
formation process of all massive galaxies. 
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